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scale and scheme. Any perturbative series in QCD is shown to be equivalent 

^ , to a particular point in this function. This function can be computed from a 
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, set of first-order differential equations involving the extended beta functions. 

We propose the use of these evolution equations instead of perturbative se- 
ries for numerical evaluation of physical observables. This formalism is free 
of scale-scheme ambiguity and allows a reliable error analysis of higher-order 
corrections. It also provides a precise definition for Aj^ as the pole in the 
associated 't Hooft scheme. A concrete application to R{e^e~ hadrons) is 
presented. 
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The scalc-schcmc ambiguity problem [1,2,3] remains as one of the major cor- 
nerstones impeding precise QCD predictions. Although all physical predictions in 
QCD should in principle be invariant under change of renormalization scale and 
scheme, in practice this invariance is only approximate due to the truncation of 
their perturbative series. 

Consider the A'^-th order expansion series of a physical observable R in terms 
of a coupling constant as{lJ>) given in scheme S and at a scale /i: 

Rn = roa^sin) + ri(/i)a^+\/i) H h rN{ij)a^^^ {jj) . (1) 

The infinite series i?oo is renormalization scale-scheme invariant. However, at any 
finite order, the scale and scheme dependencies from the coupling constant as{lj) 
and from the coefficient functions ri[ii) do not exactly cancel, which leads to a 
remnant dependence in the finite series. Different choices of scale and scheme then 
lead to different theoretical predictions. The availabihty of next-to-next-to-leading 
order results in QCD [4,5] has accentuated the need for study on the scale-scheme 
dependence (e.g., see Ref. [6]). 

There have been traditionally two positions on this subject. The first one is 
to consider the scale-scheme ambiguity as intrinsically unavoidable, and interpret 
the numerical fluctuations coming from different scale and scheme choices as the 
error in the theoretical prediction. This point of view, aside from being overly 
pessimistic, is also very unsatisfactory. First of all, in general we do not know how 
wide a range the scale and scheme parameters should vary in order to give a correct 
error estimate. Secondly, besides the error due to scale-scheme uncertainties there 
is also the error from the omitted higher-order terms. In such an approach, it is 
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not clear whether these errors are independent or correlated. The error analysis in 
this context can become quite arbitrary and unreliable. 

A second approach is to optimize the choice of scale and scheme according to 
some sensible criteria. Commonly used scale setting strategies include the Principle 
of Minimum Sensitivity [1] (which also optimizes the choice of scheme), the Fastest 
Apparent Convergence criterion [2] and the BLM method [3]. 

In this paper we propose the use of the Extended Renormalization Group 
Equations as a transparent solution to the scale-scheme ambiguity problem. In 
this approach, a perturbative series only serves as an intermediate device for the 
identification of scale and scheme parameters. The ultimate prediction is obtained 
through evolution equations in the scale- and scheme-parameter space. This ap- 
proach sets the ground for a reliable error analysis and also provides a precise 
definition for Aj^. 

We will consider the case of QCD with Nf massless quarks. Let us first explain 
the concept of the universal coupling function in QCD. It is well-known that at the 
renormalization stage a particular subtraction prescription and a particular renor- 
malization scale must be specified. Let us parametrize the subtraction prescription 
by an infinite set of continuous "scheme parameters" {ci} and the renormalization 
scale by The universal couphng function (see Fig. 1) is the extension of an 
ordinary coupling constant to include the dependence on the scheme parameters 

{a}: 

a = a(/i/A, {a}) . (2) 

For the moment, the presence of the quantity A in the previous expression can be 
justified on dimensional grounds. We will identify it later with the 't Hooft scale 
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of the chosen scheme. 

Stevenson [1] has shown that one can identify the beta-function coefficients of 

a renormaUzation scheme as its scheme parameters. That is, if a given scheme has 
the following beta-function expansion: 

then, the coefficients {Pi, i = 2,3, . . .} can be considered as the corresponding 
scheme parameters. (The scheme invariance of the first two coefficients /3o and /3i 
is a well-known fact. It is important not to confuse the coefficients /3j with the 
functions to be introduced later.) It will be very convenient to use the first 
two coefficients of the beta functions to rescale the coupling constant and the scale 
parameter log/x^. Let us define the rescaled coupling constant and the rescaled 
scale parameter as 

Then, the rescaled beta function takes the canonical form: 

/3(a) = ^ = -0^(1 + + 020^ + 030^ H ), (5) 

or 

with Cn — PnPo~^/Pi for n = 2, 3, • • •. This rescaling process serves to "unitarize" 
the expansion coefficients. For a well-behaved scheme in QCD, we would expect 
its beta-function expansion to roughly resemble a geometrical series, at least for 
the first few coefficients. In fact, for the MS scheme we have — j3^'^j3Q/0l , 
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where [7] 

/3o = 11 - ^Nf , 

A = 102 - ^Nf , (6) 

^MS 2857 5033 ^^ 325 ^. 

2 18 ^ 54 ■'^ ' 

and for Nf = 0,1,2,3,4,5,6 we have respectively cf^ = 1.5103,1.4954,1.4692, 
1.4147, 1.2851,0.92766, -0.33654. We can clearly see that indeed cf^ is of order 
of magnitude unity. This should be contrasted with the large value of /^l^^ before 
the rescaling process. 

The Extended Renormalization Group Equations simply express the invariance 
of physical quantities under scale- and scheme-parameter transformations. These 
equations have been studied long ago by Stiickelberg and Peterman [8] and later 
by Stevenson [1]: 

5R BR OR ^ 

or oa or 

_ dR dR _^ ^ ' 

5cn da dcn 

Various quantities in these equations need explanation. 

First of all, notice the distinction between 5/Sr and 8/ dr. The partial deriva- 
tive S/Sr takes into account the full variation of R under r transformation, whereas 
the partial derivative d/dr affects only the expansion coefficients of R. An analo- 
gous distinction holds between 5/5cn and d/dcn- In other words, in the left-hand 
side 

R = i?(T, {a}) , (8) 

whereas in the right-hand side 

R^R{a,rn{r,{ci})) . (9) 



Mathematically these are two different functions since their domains are different. 
However, it is common practice to use the same notation for both and draw the 
difference only for their partial derivatives. 

The fundamental beta function that appears in Eqs. (7) is defined as: 



Set 

/3(a, {ci}) = 7^ = -0^(1 + a + c^a^ + 030^ H ) , (10) 

or 



and the extended or scheme-parameter beta functions are defined as: 



/3(„)(a,{c,})^-^. (11) 



As shown by Stevenson [1], these extended beta functions can be expressed in 
terms of the fundamental beta function. Indeed, the commutativity of second 
partial derivatives 



implies 



SrScn ScnSr 



St 5cn 



(12) 



(13) 



where /?^^^ = 5/?(„-)/9a and P' — 8/3 /da. Prom here 



PPL) = P{n)P' - , (14) 



2 / Pin)\ _ n+2 



%)(«> {ci}) ^ -P{a,{ci}) J dx 



a 

j.n+2 



(15) 







where the lower hmit of the integral has been set to satisfy the boundary condition 



That is, a change in the scheme parameter can only affect terms of order a"'"^ 
and higher in the evolution of the universal coupling function [1]. 

The extended renormalization group equations in (7) can be put into the form: 



These equations can now be interpreted in the following manner. The left-hand 
sides represent the variation of the expansion coefficients of R under scale-scheme 
transformations. Thus, a given perturbative series can be evolved into another 
perturbative series, provided we know the extended beta functions on the right- 
hand side. In fact, in the expansion series of it! as given in Eq. (1), the only 
quantities that cannot be modified under scale-scheme transformations are the 
tree-level coefficient tq and exponent p, since they are renormalization scale- and 
scheme-independent . 

We can standardize a perturbative series R = roa^ + - • ■ by defining an effective 
charge aji (see Grunberg in Ref. [2]) by: 



Since R, ro and p are all renormalization scale and scheme invariant, the effective 
charge is also scale and scheme invariant. We only need to study the evolution 
of one effective charge to another. The appropriate values of ro and p can always 
be put back at the end of the analysis. 




(16) 




dcn ^("^aa ■ 



(17) 




(18) 
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More concretely, given two perturbative series: 

R = roa^ = roaP + ria^+^ H , 

(19) 

R' = r'oa^, = r'oaP' + r'iaP'+^ + ■ ■ • , 
we can evolve R into R' by following the next three steps 

1) Obtain the effective charge of R 

a^^(^Y' = a + ^a^+(^+'-^ 'i)a> + ... . (20) 
\roJ pro \pro 2 r^J 

2) By using the extended renormalization group equations, we can change the 
perturbative coefficients of aR and evolve it into 

= a + + + ■ (21) 

p' r'o \p' r'o 2 p'^ r'l ) ^ ' 

3) Obtain R' from aR' by 

R! = tWr, . (22) 

Hence, by using the extended renormalization group equations and by appro- 
priately raising the power and rescaling our final result, we can evolve a given 
perturbative series R into any other perturbative series R'. The commutativity 
of second partial derivatives guarantees that the final result is independent of the 
path chosen for evolution. 

Since the universal coupling functions a(r, {cj}) covers all possible choices of 
scale and scheme, any effective charge qr can be expressed in terms of it: 

aR = a{TR,{cf}) , (23) 

where tr and cf' are respectively the scale and scheme parameters of R. Notice 
that this is also true for multiple-scale processes: given the perturbative series of a 
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multiple-scale process, we can also define an effective charge associated to it, and 
this effective charge can then be written in terms of the universal coupling function 
at a particular set of values of r and {cf^}. 

The universal coupling function is dictated by the evolution equations: 

^ = /3{a, {ci}) = -a^(l + a + c^a^ H ) , 

or 

Sa } (2^) 

^ = %)(a, {c^}) = -Pia, M) J dx^^^^^ . 



We shall define a(T, {cj}) here with the boundary condition: 

a(0,{0}) = oo. (25) 

Notice that the above equations contain no explicit reference to QCD parame- 
ters such as the numbers of colors or the number of flavors. Therefore, aside from its 
inflnite dimensional character, a(r, {q}) is just a mathematical function like, say, 
Bessel functions or any other special function. Truncation of the fundamental beta 
function simply corresponds to evaluating a(T, {c^}) in a subspace where higher 
order Ci are zero. In principle, this function can be computed to arbitrary degree 
of precision, limited only by the truncation of the fundamental beta function.* 

Notice that the 't Hooft scheme [9] deflned by a'tui'T) = ^(t, {0}) is privileged 
since it is totally devoid of higher-order corrections. In fact, a'tH(T) is exactly given 
by the solution of 

1 +log(^^)=r. (26) 



Q'tH V 1 + O'tH 

For any single-scale process R{ij) there exists a scale n — for which the 



* The n\ divergence expected for the beta function coefficients of physical effective charges 
can impose some theoretical limit to the achievable precision. Further discussion in this 
direction is beyond the scope of this paper. 



scale parameter tr = 2/3q/3{'^ log(/x/A^^) vanishes. We will call A^^ the 't Hooft 
scale ^ of the i?-scheme. To understand the meaning of the 't Hooft scale, let us 
consider the MS scheme coupling constant: 

«Ms(/^) - « log(/^/AS)' if"}) ■ (27) 

Notice that a priori we do not know the behavior of a^{fi) at /i = A^^: it could 
be infinite, finite, or simply not well-defined. However, fj, = is the pole in the 
't Hooft scheme associated* to the MS scheme: 

«'tH-Ms(/^) ^ « log(/^/Am)' {0>) ' (28) 

because a(0, {0}) = oo by boundary condition. Since the 't Hooft scheme is com- 
pletely free of higher-order corrections, this provides a precise definition for Aj^. 
This cures the well-known arbitrariness in the definition of Aj^ due to the presence 
of higher order corrections in the MS scheme. 

Given two effective charges, say 

«ms(/^) = «(^M/^MS)'{^F^>) 
aRiQ)^a(^^logiQ/Al''),{c^}^ , 

we can expand aR{Q) in a power series of Oy^^ii). To this end, we can perform a 

t For multiple-scale processes, the sub-manifold where the scale parameter vanishes defines 

the '"t Hooft surface" . 
★ There are infinite 't Hooft schemes, differing only by the value of the 't Hooft scale A 

The word "associated" here means we are choosing the particular 't Hooft scheme that 

shares the same 't Hooft scale with the MS scheme: A *^ = A^. 

MS 
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(29) 



Taylor expansion around the point 



Pi 



(30) 



and generate a series representation for aji{Q): 



adQ) = a[^ log(/x/AS) + f, {cf^ + Ci} 



Pi 



MS 

6a 



+ 
+ 



1 

2! 
1 
3! 



r6^\ _2 ^ 2 

ySrScn 



TCn + 



SCnSCr 



(31) 



where 



T 



6^a 



2JI 
Pi 



M MS 



+ 



(32) 



, 



and a and its derivatives are evaluated at the point specified in Eq. (30). To order 
a^, we only need the following partial derivatives: 



/ 5a 

5a 

5C2 

/ 5a 

5^ 
5r'2 
_fa_ 
5t5c2 
5^ 
5r^ 



P^-a^ -a^- C2a^ + O(a^) , 

/3(2) = a3 + 0(a5) , 

%) = ^«' + 0(a5), 

2a^ + 5a^ + 0(0^) , 

-3a^ + O(a^) , 

-6a^ + O(a^) . 



(33) 
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After grouping all the terms in powers of a = we obtain: 

O'RiQ) = a — fci^ + (c2 — f + f ^) 

+ Qc3 - (C2 + 3c2)f + - f 3^ (34) 
+ 0{a^) , 

where f and Cn are as given in Eq. (32), and C2 = cf^. Although we have used R 
and MS schemes in the derivation of this last formula, naturally it is also valid for 
any other pair of effective charges. Notice the occurrence of f and q in all higher 
order coefficients. By using the evolution equations in (24), we are effectively 
performing a partial resummation of the perturbative series to all orders. 

Any physical quantity R calculated perturbatively in the MS scheme can be 
put into the above standard form. Hence, we can use Eq. (34) to identify the scale 
and scheme parameters of the effective charge of R. In other words, to give the 
numerical prediction for R — r^aP + • • • we proceed to 

1) compute its perturbative series in terms of the coupling constant in some 
scheme, say MS scheme, 

2) standardize the series and identify the scale and scheme parameters {tr, {cf}) 
order- by-order via Eq. (34), 

3) from the knowledge of evolve the universal coupling function a(r, {cj}) 
to the point [tr, {cf}). Put back the values of tq and p if necessary. 

Naturally we can also go in the opposite direction and obtain A^^ from the 
experimental measurement of R. In Fig. 2 we show the various experimental and 
theoretical errors involved in the analysis. For the measurement of A^, the input 
experimental error must be combined with the scheme uncertainty to give the error 
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estimate for A^^. Similarly, for the prediction of aji^r) the error from must 
be combined with the scheme uncertainty in order to give the prediction error. 

As a concrete example, let us consider the total hadronic cross section in e'^e~ 
annihilation R{Q) = R{e^e~ — > hadrons) recently calculated to order [4,5]. For 
five light-quark flavors we have 



R{Q) 



11 

y 
11 



1 + - + 1.4092 

TT 



12.8046 



1 + 



TT 



where we have used a = a^^^Q) for the strong coupling constant. 
By putting R{Q) into the standard form of Eq. (34): 



(35) 



aniQ) = a+ 1.1176 - 8.05426 , 



(36) 



where a = aTju{Q), we can identify 



f = -1.1176 , C2 = -10.4209 . 



(37) 



Knowing that f = 2/3^/3f ^ log(Ajg/A^H) g^^^j = cf - cf^, we conclude that 



= 1.4443 A^ 



'tH 



MS 



= -9.4932 . 



(38) 



Experimentally [10] we have 



r(31.6 GeV) = — i?(31.6 GeV) = 1.0527 ± 0.0050 , 



(39) 



which gives 



ai?(31.6 GeV) = 0.0665 ± 0.0063 . 



(40) 



We have to take into account the scheme uncertainty in addition to the ex- 



perimental error in order to quote a correct error estimate for A^ (see Fig. 3). 



MS 
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The scheme uncertainty of R can be quantified by a reasonable estimate of its next 
scheme parameter: c^. G. B. West [11] has put an estimate r4 = —158.6 for the 
coefficient r4 in 



(41) 



By reducing this equation into the standard form of Eq. (34), we arrive to a value 
C3 — — Cg^^ = —99.474. Assuming is of order unity, we conclude that 
Icg^l ~ 100. Although we have some reservation on West's estimate (see Ref. [12]), 
we shall nonetheless use it to illustrate our procedure. A better estimate of c§ will 
lead to a better error estimate for A^j^. 

MS 

In Fig. 3 we show the universal charge for aoir) = a{T, {c2 = c^^, cs = 04 — 
... = 0}) and its evolution under a scheme uncertainty C3 = ±100 to a±{T) — 
air, {c2 = cf , C3 = ±100, C4 = C5 = . . . = 0}). 

The evolution in the scheme parameters are dictated by: 



5a 

SC2 



/4 
dx^ ~ + 0{a^) , 

(42) 

Sa f 1 



For our region of interest (a ~ 0.07) the first term in each expansion series suffices. 
But we should use the full integro-differential equation whenever we want to evolve 
a to a higher- value region. This will not only improve the accuracy of our result, 
but also will respect the commutativity of the second-order partial derivatives of a 
and thus ensure the independence of the result on the choice of integration path. 

To obtain ao and a± defined above, we follow the next steps: 
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1. — Generate the 't Hooft scheme couphng O'tH by solving iteratively 

T + log (1 + l/o/tH) 

2. — Evolve a/tH to ao by displacing in C2. Prom (42) we have: 

ao = "-^ ^2 • (44) 

(l-2cf a^tH^)'/' 

3. — Evolve ao to a± by displacing in C3. From (42) we have: 

a± = TTs • (45) 

In Fig. 3 we show the various errors involved in this analysis. Numerically we 
find the experimental, scheme, and total errors for r to be: 

Arexp = (t6 - T2)/2 - 1.41 , 

Arsch = (rs - T3)/2 = 0.22 , (46) 
Artot = (ty - Ti)/2 = 1.63 . 
These errors can be translated into uncertainties in since there is a one-to-one 

MS 

correspondence between r and We can see that most error comes from the 

experimental error in aji. We can also see that the experimental error and the 
scheme error are highly uncorrelated since Artot ~ Arexp + Argg/j. Numerically we 
obtain ri = 10.129, r4 = 11.666 and — 13.379. Knowing that 

2/32 / 31.6 GcV \ 

we arrive at the following result for A^: 

° MS 

= 470igg MeV . (48) 
If there were no experimental error, the estimated scheme uncertainty would lead 
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toAg = 470l^S MeV* 

As a second application of our formalism, we will show next how to use the 
experimental result of 07^(81.6 GeV) to predict other effective charges. Specifically, 
we will give a prediction for a^{Mz), where Mz = 91.173 GeV is the mass of 
the Z-boson. The evolution of a_R(31.6 GeV) to a^{Mz) is illustrated in Fig. 
4. Notice that the experimental and the scheme uncertainties confine the correct 
result for qr into an approximate parallelogram ABCD. We can then evolve 
this parallelogram into any other scheme and scale. We will use — ±100 and 
c|^S _ _|_-|^ estimate the scheme uncertainties in aR and Oj^. For aj^(M^), 
the parallelogram ABCD is evolved into the parallelogram A'B'C'D' . Notice the 
inversion of the orientation of the new parallelogram due to the opposite signs of 
C2 and (^^. Notice also the absence of scheme uncertainty in the 't Hooft scheme. 

From TA — 10.129 and tq — 13.379 and knowing that — tr — f with 
f = 2/3^/3f 1 [log(M^/Aj^) - log(31.6 GeV/A^j)] = 4.339, we find ta' = 14.468 
and Tc = 17.718. From here and using = 0.92766 and c^^^ = 1 in the MS 
version of the Eqs. (44) and (45), we obtain: a^' = 0.05772 and a'^' = 0.04818. 
Hence, we arrive at the prediction 

a^{Mz) = 0.0530 ± 0.0048 , (49) 

or equivalently,^ 

-k We have roundc^d off the above results. More precise values are = 472^204 MeV 

for the case including experimental error, and = 4721:11 MeV for the case without 
experimental error. 

t This value is higher than the world average aj^{Mz) = 0.1134±0.0035 quoted in the Review 
of Particle Properties [13] but still consistent with other quoted values for ay^{Mz). For 
instance, a-^{Mz) = 0.118 ± 0.008 is obtained by OPAL [14]. The detailed analysis of 
consistency between the various experimental results is beyond the purpose of this paper. 
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ajj^{Mz) = 0.132 ± 0.012 . (50) 

Let us show next that our formahsm is closely related to the FAC (Fastest 
Apparent Convergence) criterion when only the next-to-leading-order coefficient is 
known. We define FAC here as the condition of a vanishing next-to-leading-order 
coefficient. 

Given 



with 



f = 2PiP^^ [log(g/A^H) - log(/x/Ag)l , (52) 



and assuming a complete lack of knowledge of the scheme parameters {c^^} and 
{cf}, we cannot do much better than to approximate: 



MQ) = «(^iog(Q/A^"),{0}) , 



(53) 



2j^0 /A'tH 



and absorb the uncertainty from scheme-parameters into our theoretical error. 
However, the last two equations imply that: 

/ AM \ 

aR{Q) = a^U=j^Q\ . (54) 

Hence, by setting the coefficient f in Eq. (51) to zero, we obtain the correct ratio 
for the two scales. Actually, this holds true in general: by applying the FAC 
criterion, we always obtain the correct ratio of the 't Hooft scales, whether or 
not we know the scheme parameters. Thus, despite its apparent naiveness, FAC 
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actually constitutes a correct first step towards the elimination of scale-scheme 
ambiguity. 

Finally let us comment on the definition of effective charges. The definition of 
aR hy R = roa^ is not the only possibility [15]. In fact, we can apply the extended 
renormahzation group technique to the effective BLM charge [3] in the following 
manner. Given two physical quantities R and R' computed in a particular scheme 

R{Q) = ro{Q)aP{ii) + (rio(g) + rn(g, /.)/3o) aP+\i^) + • • • , 

(55) 

R'iQ') = r'M)aP\^i') + (rlo(g') + r'n(Q', ^i')Po) of>'+\^^') + • • • , 

where (Sq = 11 — |A^j is the first beta function coefficient, we can define their 
"effective BLM charges" by 

R{Q) = ro(Q)c4_5i^(g) + rio(g)c4+ slm(Q) + " ' ' > 

(56) 

R'iQ') = r'^{Q')oFR>_BLMiQ') + r'M')4?'-BLMiQ') + • • • . 

with olr-blm{Q) — OiifJ'*) where /i* is the solution of rii{Q, //*) = 0, and similarly 
<^R'-BLm{Q') — ck(a*'*) where is the solution of r[i{Q' , /j,'*) — 0. With this 
choice of scale, vacuum polarization contributions are associated with the charge 
rather than the expansion coefficients, and the scale tends to reflect the mass of 
the virtual gluons. We can then apply the evolution equations to o-R-BLuiQ) 
and evolve it to olri_blm{Q')- Alternative definitions of effective charges have 
also been discussed in Ref. [15]. However, as we have shown in this paper any 
convenient choice of effective charge can be used to relate physical observables. In 
practice we can adhere to the definition R = roa^ which has the advantage that 
ro and p are renormahzation scale-scheme invariant quantities. 
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To summarize, we have explained the use of extended renormahzation group 
equations to relate physical observables. The most distinctive feature of this for- 
malism is that, in this approach, the perturbative series of a physical observable 
only serves to identify the scale and scheme parameters. The final prediction is 
obtained by the evolution of a universal coupling function. The prediction is scale- 
scheme independent in the sense that given the initial perturbative series in any 
scheme at any scale, we will always obtain its correct scale and scheme parameter 
and hence arrive at the same prediction. We have shown that this formalism sets 
the ground for a reliable error analysis, and that can be unambiguously de- 
fined as the pole in the associated 't Hooft scheme. Finally, we have shown that 
this formalism is equivalent to the fastest apparent convergence criterion in the 
absence of information on scheme parameters. 
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FIGURE CAPTIONS 



1) Pictorial representation of the universal coupling function a(r, {q}), where 
T is the scale parameter and {cj} the scheme parameters (see definitions in 
text). 

2) Graphical representation of the various errors involved. For the measurement 
of (or equivalently A^^) the input experimental error must be combined 
with the scheme uncertainty. For the prediction of aji{T), the error in 

IVio 

must be combined with the scheme uncertainty. 

3) Measurement of A^ from the experimental result of ajj(31.6 GeV) . We have 
parametrized the scheme uncertainty with a value c§ = 100. The scheme, 
experimental and total errors are respectively given by Agch — ("^5 ~ ''3)/2, 
^exp = {'T6 — T2)/2 and Atot = ("77 — ''■i)/2. There is a one-to-one relationship 
between r and aJ^ given by r = 2/?^/?rMog(31.6 GeV/1.4443 A^). 

4) Prediction of ajj^g{Mz) from the experimental result of aij(31.6 GeV). By 
using the extended renormalization group equations, the quasi-parallelogram 
ABCD is evolved into the quasi-parallelogram A'B'C'D' . Notice the inver- 
sion of the orientation of the parallelograms due to the opposite signs of 

and Cg'^. Notice also the absence of scheme uncertainty in the 't Hooft 
scheme. 
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